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The ATLAS and CMS experiments have recently claimed discovery of a Higgs boson- 
like particle at <~ 5cr confidence and are beginning to test the Standard Model predictions 
for its production and decay. In a variety of supersymmetric models, a neutralino NLSP 
can decay dominantly to the Higgs and the LSP. In natural SUSY models, a light third 
generation squark decaying through this chain can lead to large excess Higgs production 
while evading existing BSM searches. Such models can be observed at the 8 TeV LHC in 
channels exploiting the rare diphoton decays of the Higgs produced in the cascade decay. 
Identifying a diphoton resonance in association with missing energy, a lepton, or 6-tagged 
jets is a promising search strategy for discovery of these models, and would immediately 
signal new physics involving production of a Higgs boson. We also discuss the possibility 
that excess Higgs production in these SUSY decays can be responsible for enhancements of 
up to 50% over the SM prediction for the observed rate in the existing inclusive diphoton 
searches, a scenario which would likely by the end of the 8 TeV run be accompanied by 
excesses in the 77 + £/MET and SUSY multi-lepton/& searches and a potential discovery in 
a 77 + 2b search. 
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I. INTRODUCTION 

The ATLAS and CMS experiments at the Large Hadron Collider (LHC) have recently reported 
discovery of a new particle consistent with the Higgs boson of the Standard Model [HE]. Each ex- 
periment independently reports nearly 5<t combined local significance for a mass of rrih ~ 125 GeV, 
largely fueled by the searches for resonances in diphoton and ZZ* —> 4/ final states. 

The discovery of the Higgs marks the final chapter in the quest to discover the particle content 
of the Standard Model. It also represents the start of a new era of Higgs physics at the LHC. With 
additional data and further analyses the experiments will be able to study the properties of this 
new particle, in particular its production modes and decay branching ratios. The Standard Model 
makes very definite predictions for these properties (up to calculational uncertainties) given the 
mass of the Higgs, so any significant deviation from the expected values would be an immediate 
sign of new physics. 

There are many ways that physics beyond the Standard Model (BSM) can manifest itself in 
Higgs events. One possibility is that the Higgs couplings to known SM particles are modified by new 
physics, such as loop effects of new particles, mixing of the Higgs with other new states (extended 
Higgs sectors), and Higgs compositeness. Such effects will modify both the decay branching ratios 
and production modes of the Higgs, potentially giving observable deviations from SM predictions. 
Several works perform general analyses of these possibilities in the context of recent LHC Higgs 
results [Ml]. 

A second possibility is that on-shell production of BSM particles can lead to Higgs production 
in the decay cascade. This possibility has been discussed as a discovery mode for an SM-like Higgs 
boson in several models |12H20j . for the extra Higgs states of the minimal supersymmetric standard 
model (MSSM) and its extensions [2TH28] . and for composite Higgs models [29ti32] . Conversely, the 
discovery sensitivity for BSM physics through Higgs production has been studied in the context 
of the MSSM with a bino LSP [33H35] . a gravitino LSP [36H4"0] . and additional feebly coupled 
LSPs |41H43| . as well as in the context of vector-like heavy quarks |3 HI44fH5] . Higgs production in 
cascade decays has also been studied as a way to probe the details of the the spectra and couplings 
of BSM particles after their discovery |49H54j . 

In this work we consider a class of SUSY models with neutralino NLSPs which dominantly decay 
to a Higgs boson. In particular we focus on models with light third generation squarks, which are 
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motivated by naturalness of the electroweak scale. Such models can remain unconstrained by 
existing SUSY searches while still having total SUSY production cross sections at the 7 or 8 TeV 
LHC of up to a few picobarns- comparable to the total SM Higgs production cross section and 
significantly larger than the cross section for associated production in the SM. The sensitivity of 
SUSY searches is limited due to the suppression of missing energy when a light NLSP (< 200 GeV) 
decays to a Higgs and the LSP, as has been explored previously in more general models [67 \ 142 1 14*3"]. 

The dominant decays for a light Higgs are h —> bb and h — > WW*, and much of the existing 
work on BSM Higgs production has focused on these modes. The h — )• bb mode suggests search 
channels with many 6's or searches for 'Higgs-tagged' boosted jets [El ESJ EZl EQj [35], while the 
h — > WW* decay suggests searches in leptonic channels with missing energy. The many-6 modes 
have large QCD backgrounds when the MET signature is suppressed, while the leptonic modes 
have smaller branching fractions and can have significant backgrounds as well. Furthermore, it can 
be difficult to link a new physics signal with the Higgs in these channels, especially in the presence 
of other hard objects from the SUSY cascade. 

We will focus instead on the potential for observing such models using the rare but distinctive 
diphoton decay mode of the Higgs, which has a branching fraction of 0.22% for a 125 GeV SM 
Higgs. We show that searches requiring a diphoton with invariant mass near 125 GeV plus other 
objects such as 6-tagged jets, missing energy or leptons can target the models we consider while 
having very low SM backgrounds (a similar point is made in the context of searching for heavy 
vector-like quarks in ref. |31j). Since such searches look for a narrow resonance in the diphoton 
mass, the backgrounds can be directly determined from data and the observation of a peak near 
125 GeV would strongly suggest new physics involving production of an on-shell Higgs. We find 
that in the 8 TeV LHC run the 77 + 2b channel could give a strong signal in regions of parameter 
space where other searches are not sensitive, and that the 77 + MET and 77 + I channels have 
sensitivity comparable to and in some regions exceeding that of existing SUSY searches. These 
results are summarized in fig. [6j 

Although targeted 77 + X searches have the best potential to probe our models, for sufficiently 
light squarks even inclusive searches for the SM Higgs boson in the resonant 77 and ZZ* — > 41 
channels can have some sensitivity. Enhancements of 20-50% over the SM rates can be accom- 
modated within current bounds from other searches. Searches in non-resonant channels, such as 
h — > WW*, h — y tt and Vh — > Vbb, generally include additional cuts that tend to reject SUSY 
events for the models we consider, so no excess rate is predicted for these analyses. Our models 
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therefore predict a very particular pattern of excesses in certain Higgs searches which can be ob- 
servable with sufficient data. Such excesses can be consistent with the currently observed rates in 
the ATLAS and CMS Higgs searches, which report a higher than expected rate in the diphoton 
channels. We find that in the regions of parameter space where the excess rate in SM Higgs searches 
is significant, strong signals are expected in the 77 + 6/MET modes within the 8 TeV run. 

We now give an outline of the rest of this work. In sec. [H] we describe models in which a 
neutralino NLSP can decay dominantly to the Higgs and in which the third generation squarks 



are light and give the dominant SUSY production cross section. In sec. Ill we consider the limits 
on these models from existing LHC SUSY searches. We then consider the signal in the Higgs 



to diphoton channels (sec. IV), determining the sensitivity of existing searches for 77 + i and 
77 + MET, and suggesting a new resonant search for diphotons in association with a 6-tagged 
jet. The possibility that current hints of excesses in the inclusive Higgs diphoton channel can be 
explained within this model is discussed in sec. [Vl We also briefly discuss the signals from direct 



electroweak production of neutralinos (sec. VI). We then conclude in sec. VII with an outlook for 
the rest of the 8 TeV run and a discussion of more general motivations for searches for new physics 
in the Higgs + X channel. 



II. MODELS 

In this section, we describe a set of models in which neutralino NLSPs decay dominantly to the 
LSP and a Higgs. An important consequence of this mechanism for Higgs production is that if the 
LSP is light and the splitting between the NLSP and the Higgs is not too large, the MET signature 
in the SUSY cascade is suppressed. This can relax the collider constraints on superpartners, as 



we will discuss in detail for our model in sec. Ill and is explored for gluino and degenerate squark 



simplified models and the c(N)MSSM in refs. |37} W2\ I43j . In particular, this can allow for large 
colored production cross sections with only weak constraints from existing searches. 

Amongst the SUSY particles the squarks and gluinos typically have the largest production cross 
sections due to their color charge. If the squarks are degenerate in mass, existing SUSY searches 
strongly constrain their mass in the simplest .R-parity conserving scenarios. Squarks decaying 
directly to light neutralinos must be heavier than ~TeV |55|, I56j . These bounds are in tension 
with naturalness of the electroweak scale in SUSY, as if the stops are heavy then they give large 
corrections to the Higgs masses that must be tuned away. 
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It is possible however to recover SUSY naturalness if only the third generation squarks are 
light, with masses ~ few hundred GeV, with the other squarks much heavier- so-called "natural 
spectra" |57H7U| . (Some extension of the MSSM is then necessary if the lightest Higgs is to have a 
mass near 125 GeV; a number of possibilities have been explored |71H88j .) This scenario is much 
less constrained than the case of degenerate squarks since the third-generation squarks have smaller 
production cross sections and can produce less missing energy in their decays 



Combining these two elements, a natural squark spectrum and a neutralino NLSP decaying 
dominantly to the Higgs, we find a set of models that evade constraints from existing searches 
but give large cross sections for Higgs production, potentially giving striking signatures in searches 
exploiting the Higgs resonance. The key features of these realistic models can be cast into a 



simplified model (section II C ) which will be our main focus for the rest of this work. 



A. NLSP Decays to the Higgs 

If imparity is conserved, the two-body decays of a neutralino NLSP to a neutral LSP can produce 
an on- or off-shell Higgs, 2f, or photon. In this section we describe some types of models in which 
a wino, bino, or higgsino NLSP decays dominantly to an on-shell Higgs and a neutral LSP, making 
this channel a promising way to probe such a model at the LHC. 

States uncharged under the Z 

At the renormalizable level, one simple scenario for obtaining large branching ratios of the 
NLSP to Higgses is to have both the LSP and NLSP be formed mostly from states that do not 
couple to the Z or the photon (i.e. states with T 3 = Y = 0). This is achieved in the MSSM with 
a neutral wino-like NLSP and a bino-like LSP (as studied in |33H35| ). When \x 3> M\,M.2, the 
mixing of the NLSP and LSP with the higgsino is small. Decays of the NLSP to the Higgs are 
suppressed by one power of this mixing, while decays to the Z require mixing of both the NLSP and 
LSP states and are suppressed by two powers of the mixing. This is illustrated diagramatically in 
fig. [Tj From an effective operator viewpoint, integrating out the higgsinos generates the dimension 
five operator (h^a 3 h)W°B, while decays to the Z will arise only from the dimension six operator 
(since the dimension five operator does not involve couplings to the neutral 
goldstone mode of the Higgs) and are therefore suppressed by additional powers of the higgsino 
mass. The dimension five operator Wa^ u BW^ u can be generated by loops involving sfermions and 



allows decays to the photon or Z, but the loop factor suppression renders this negligible. Fig. 2(a) 
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FIG. 1: Diagrams for decays of a wino NLSP to a bino LSP and a Higgs or a Z. The neutralinos are shown 
in the gauge eigenbasis, with mixing represented by interaction vertices (crosses). The decay of the wino to 
the bino and a Higgs proceeds through mixing of one of the states with the Higgsino, while the decay to the 
Z requires at least two mixings and is further suppressed. 

shows the branching ratio of a 175 GeV wino to a Higgs and a nearly massless bino as a function of 
[A and tan/3, for a Higgs mass of 125 GeV in the decoupling limit. Though large tan/3 suppresses 
the wino decay to the Higgs, it can be highly dominant even for relatively light higgsinos. 

In extensions of the MSSM there can be additional states uncharged under the Z but with 
renormalizable couplings to the Higgs and higgsinos. This can cause the NLSP to decay dominantly 
to the Higgs by the same mechanism as above. Examples include the singlino S of the NMSSM |42j 
and string photini E2]. A singlino LSP mixes through the operator XS(H u Hd + HdH u ), and 
so unlike the bino has a decay width to the Higgs that increases with large tan /3 and is suppressed 
by low tan/3. String photini are superpartners of light U(l) gauge fields from string theory [92], 
which can have small mixings with the bino. A string photino LSP therefore couples to the Higgs 
exactly as the bino does but with a coupling suppressed by the wino-photino mixing e. A similar 
enhancement of the Higgs branching ratio can also occur for decays of a wino or bino NLSP to the 
goldstino of an it!-symmetric SUSY breaking sector |41j . 

Since the LSP can be weakly coupled to the MSSM fields in these scenarios, we will consider the 
limit in which all of the other MSSM fields decay promptly to the bino or neutral/charged winos, 
which then decay to the LSP. We restrict our attention to the case where the decay to the LSP is 
not displaced on detector scales (cr < 0.1 mm), ensuring accurate reconstruction of collider events 
and avoiding constraints from searches for metastable particles. 

Higgsino NLSP and gravitino LSP 

Another scenario for obtaining large NLSP branching ratios to the Higgs can occur for a 
higgsino- like NLSP decaying promptly to a gravitino LSP in models with low-scale SUSY breaking. 
Refs. [37H4Tj] discuss this possibility and some of the collider signatures at the LHC and Tevatron. 
If sign(^(|| 1 + tan 2 9w)) < 0, then as tan /3 approaches unity the coupling of the lightest higgsino 
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FIG. 2: (a): Branching ratio of a neutral wino NLSP to a Higgs and a bino LSP as a function of fj, and 
tan/3, for a wino mass of 175 GeV and a Higgs of mass 125 GeV in the decoupling limit. The branching 
ratio is computed using tree-level neutralino masses and mixings, with the input M\ = and choosing M% 
such that the wino mass is 175 GeV. The resulting bino mass varies throughout this parameter space but 
remains very light (< 10 GeV). (b): Branching ratio of a higgsino NLSP to a Higgs and a bino LSP as 
a function of higgsino mass and tan/3, for large Mi, M 2 , a light gravitino, and a Higgs of mass 125 GeV 
in the decoupling limit, using the formulae of [33]. This plot assumes sign(/i(^ + tan 2 9w)) < 0, which 
causes the higgsino- Z-gravitino coupling to vanish at tan /3 = 1. In both plots, the gray dashed line indicates 
tan /3 = 1.4, the approximate lower bound for perturbativity of the top Yukawa coupling at the unification 
scale ~ 10 16 GeV. 

to the Z approaches zero [HU]) due to cancellation of the contributions from the H u and com- 
ponents. The higgsino LSP can then decay dominantly to a Higgs and gravitino, as illustrated in 
fig. |2(b)| Since the gravitino LSP can be very weakly coupled, we will again consider the limit when 
all other MSSM fields decay promptly to the NLSP, which then decays promptly to the gravitino 
(requiring <J(F) < 100 TeV). 

Spectra 

In the models discussed above the NLSP is often accompanied by other nearly degenerate states 
that are important for phenomenology. We identify three classes of spectra (illustrated in fig. [3]) 
that capture the possible signals: 
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FIG. 3: Three types of spectra for the NLSP, LSP and associated states such that the NLSP decays 



dominantly the Higgs. See the text (sec. II A) for discussion of each case 



(a) Bino NLSP If the bino is the NLSP, then it can decay dominantly to a Higgs and singlino (or 
string photino etc.) LSP provided that [i is large. If the LSP is very weakly coupled, essentially 
all SUSY cascade decays will proceed through the NLSP and produce a Higgs. A pure bino 
has vanishing pair-production cross section, so this spectrum alone is not observable. 

(b) Wino co-NLSPs If the NLSP is a neutral wino with small mixings (large [i and M\ not 
too close to M2), then the charged winos are nearly degenerate with the NLSP. The decay of 
the charginos to the NLSP is then extremely suppressed; if this is the dominant decay mode 
then the charginos have a macroscopic decay length. We instead consider the case where the 
chargino decays promptly to the LSP, e.g. — > BW^ . In this scenario not all SUSY decays 
will proceed through the true NLSP (the neutral wino) and produce a Higgs. 

(c) Quasi-degenerate charginos Two types of models can give the spectrum and decays shown 
in fig. 4(c)| One is the Higgsino NLSP scenario described in sec. II A For small mixing, 
the neutral Higgsino NLSP is accompanied by a chargino and a slightly heavier neutralino. 
We will consider the minimal case where the splitting between these states and the NLSP is 
~ 1 — 10 GeV, so that they decay promptly to NLSP but the extra objects produce in the 
decay are soft and do not contribute significantly to collider signals. For a Higgsino-like NLSP 
at < 200 GeV, this is achieved with M\,Mi > 500 GeV. Because the gravitino is very weakly 
coupled, the other higgsinos will always decay to the NLSP. 

A second source for this spectrum are models where M\ « M2 so that the wino and bino are 
significantly mixed, while the LSP is an additional state such as a singlino or string photino. 
The spectrum then includes a chargino and heavier wino-bino state in addition to the mixed 
wino-bino NLSP. The splitting between these heavier states and the NLSP can be large enough 
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FIG. 4: The decay modes of light stops and sbottoms for the three types of NLSP spectra in fig. [3| See the 



text (sec. II B ) for discussion. 



to allow prompt decays to the latter (unlike the pure wino case). While the wino-bino mixing 
is suppressed at the larger values of fi required for the Higgs decay to dominate, sufficient 
mixing can still be obtained for quasi-degenerate gauginos with \M\ — Mi\ < 50 GeV. 



B. Production through third generation squarks 

SUSY naturalness motivates light stops t\ and £27 possibly with large mixing, which are generally 
accompanied by a light left-handed sbottom. When m~ 3L is small and the stop and sbottom mixing 
are negligible, the lightest squark is b^, and the tree- level splitting with ti, is given by 

ml L - ml « m| + (1 - sin 2 8 W ) cos(2/3)m| (1) 

When m^ L ~ 300 GeV, this gives a splitting mt L — mb L ~ 40 GeV. A light right-handed sbottom 
6r is not directly motivated by naturalness, but can result from UV realizations of natural spectra. 

The decay cascades of the lightest third generation squarks to a neutralino NLSP involve the 
production of a t quark (on or off shell) or b quark. When the NLSP decay produces a Higgs 
dominantly, the final states of interest will be 2h + 2(b/t) + MET. We find that final states with 
Higgses and tops are substantially more constrained by existing leptonic SUSY searches, so we 
will focus on spectra where the decays of the lightest squark(s) proceed dominantly to the 26 final 
state. 
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For each of the three NLSP spectra described in sec. we discuss the signals from the various 
third-generation squarks: 

(a) Bino NLSP If the bino is the NLSP and the LSP is very weakly coupled, then the squarks will 
decay to the bino, producing a 6 quark and, for a stop, a W boson. Sbottoms therefore give 
26 + 2h final states, while stops give 26 + 2W + 2h. The latter case is strongly constrained by 



same-sign dilepton searches (see sec. HID), requiring the stops to be heavier than ~ 400 GeV. 
A left-handed sbottom is thus generally required to be > 350 GeV, though a light right-handed 
sbottom bn can remain to give 2b + 2h at observable levels. 

(b) Wino co-NLSPs If the winos are co-NLSPs, the squarks can decay to either the charged or 
neutral wino. This decay will be dominant if the LSP is a very weakly coupled state such as 
a string photino. If the splitting between the winos and the squarks is not too much larger 
than the top mass, then decays to bottom quarks dominate, so that sbottoms decay to the 
neutral wino and stops to the charged wino. Sbottom decays through the winos then give the 
final state 2b + 2h, while stops give 26 + 2W. The latter is constrained by direct stop searches 



(sec. HID), requiring the stops to be heavier than ~ 350 GeV. 

If the LSP is a bino, then the squarks can decay directly to the LSP and a quark. This will 
be the dominant decay for right-handed squarks, so Higgses will not be produced in this case. 
Left-handed squarks will have a subdominant branching ratio to the bino (at least ~ 5%), 
but this is enough for searches in 36 + MET final states to be sensitive, as will be shown in 
sec. IIII C[ These searches constrain the left-handed sbottom to be heavier than ~ 350 GeV. 
The bino LSP scenario therefore does not have as much potential to give observable signals 
during the 8 TeV run. 

(c) Quasi-degenerate charginos As before, we will assume that the splitting between the NLSP 
and the squarks is small enough that the squarks decay dominantly to bottoms rather than tops. 
The sbottoms then decay to neutralinos while stops decay to charginos. All of these decays 
eventually produce the NLSP, possibly with additional soft decay products, so all squarks give 
the 26 + 2h final state. 
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FIG. 5: The simplified model we consider to capture a wide range of possibilities for third generation squarks 
producing Higgses in cascade decays. Only pair production of the bottom squark b is considered. 



C. A Simplified Model 



From the above discussion, we can see that the signatures of a third generation squark decaying 
to the 26 + 2/i final state can in many cases be well described by the simplified model (illustrated in 
fig. [5]) of a sbottom decaying with 100% branching fraction to a b quark and neutral NLSP which 
decays with 100% branching fraction to the Higgs and LSP, giving the final state 2h + 26 + MET. 
Some of the above scenarios give additional final state particles that carry very little energy due to 
small mass splittings. These will generally not appear as hard objects in collider events and will 
only slightly affect the kinematics of the hard 6's and Higgses. Therefore the simplified model is 
an appropriate description of the collider signatures in these models as well. 

The pair production cross sections are essentially the same for all sbottoms and stops (up to 
small electroweak effects), so this model is a good description for either a light stop or sbottom. The 
model is described by two parameters, the sbottom/stop mass and the NLSP mass. For simplicity, 
we fix the LSP mass at 1 GeV and consider a 125 GeV lightest Higgs boson in the decoupling limit, 
where its couplings are equal to those of the SM Higgs. 

We also note that several other types of models can give enhanced Higgs production with similar 
final states, i.e. Higgses + (6-quarks) + (low MET). For example, in models with additional vector- 
like generations, decays of the vector-like quarks to a Higgs and a top or bottom also produce this 
final state. In particular, decays of a new bottom-like quark to bh give the same event kinematics 
as this simplified model with NLSP mass equal to the Higgs mass, though with different production 
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cross section. The phenomenology of top-like quarks is explored in |31| . Ref. [93j showed that for 
i?-parity violating SUSY, Higgs bosons can be produced preferentially in the decays between third 
generation squarks, e.g. t% —> t\h, giving the same signatures. Another possibility is new physics 
enhancing the rare decay t — > ch, which can be observed in top pair production events [S3]. Other 
recent work studied the possibility that light colored particles in loops could increase the rate of 
Higgs pair production pp — > hh to make its effects comparable to the SM single Higgs production 
rate [95J, which would have similar signatures in multilepton and diphoton searches. We expect 
that some of our conclusions based on this simplified model will also apply qualitatively to these 



more general models, as discussed in sec. VII 



In the following sections, we describe in detail the constraints and signals of this simplified 



model in existing SUSY searches (sec. Ill) and existing Higgs diphoton searches (sec. IV). These 
results are summarized in fig. [6| In short, we find that a sbottom heavier than ~ 260 GeV with 
an NLSP lighter than ~ 200 GeV is not currently excluded. In this region there can be strong 



signals in a proposed search for resonant 77 + 1/26 (sec. IV D), and even observable enhancements 
in the inclusive Higgs diphoton searches (sees. IV A| and |V|) . We also consider the constraints on 
deviations from the simplified model, in particular non-zero branching ratios of the NLSP to the 



Z (sec. IIIB), non-zero branching ratio for the direct decay of the squark to the LSP (sec. IIIC), 



and decays of the heavier third generation squarks (sec. HID); we find that the simplified model 
is a good description of the dominant LHC signals for a wide range of realistic models. 



III. CONSTRAINTS FROM SUSY SEARCHES 



In this section we determine the constraints on models like those discussed above from existing 



searches for supersymmetry. To do so (and to determine signal rates in Higgs searches in sec. IV), 
we perform Monte Carlo (MC) simulation of collider events to determine the efficiencies of various 
LHC searches for these models. We use MadGraph 5 [SB] to generate SUSY production events, 
Pythia 6.4 |97] to decay particles, perform parton showering and hadronization and produce 
initial and final state radiation, and PGS4 [98] to simulate object reconstruction in detectors (we 
use the Pythia-PGS package to interface these tools). We use the default PGS parameter cards for 
the ATLAS and CMS detectors, but to implement 6-tagging efficiencies we use truth-level 6-jets 
and apply the p^-dependent 6-tagging efficiencies reported by the ATLAS and CMS experiments 
for each search. We scan a grid of points in Mnlsp and Mj/g spaced by 5 to 15 GeV, and generate 
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FIG. 6: The 95% excluded regions for the simplified model of sec. II C from (a) the existing searches in 

and (b) the expected 95% exclusions in 20 fb _1 at 8 TeV for the 



5 fb 1 discussed in sec. 



Ill 



and 



IV A 



same searches, extrapolated as described in the text. The diagonal hatching indicates the region constrained 
by the CMS8 36 + MET (a T ) search [P04] . horizontal for the CMS7 multilepton search [P07] . vertical for 
the CMS8 same-sign dilepton + 2b search [TUB], and solid for the CMS8 77 + MET search QTI]. The black 
solid (dashed) contours are the expected 95% exclusions in 5 fb -1 and 20 fb^ 1 from a proposed search in 



77 in association with at least two (one) 6-tagged jets, described in sec. IV D 



50,000 events at each point for the diphoton channels and 250,000 events for the non-diphoton 
channels to obtain adequate statistics for low rate searches. We obtain the cross sections for 
production of third-generation squarks from and for neutralino/chargino direct production 
from Prospino 2.1 [100J; in both cases we use the central value of the prediction and do not 
account for the theoretical uncertainty (giving more conservative bounds). We take SM Higgs 
boson cross sections and branching ratios from [TOl]. When available we use the 95% CL S visible 
cross section limits presented in the searches we study; when these are not available we compute 
the 95% CL S limit from the observed event rate and backgrounds as described in [102 . Where not 
otherwise noted, we estimate the expected limits in 20 fb -1 by linearly extrapolating the expected 
background rate and systematic uncertainty from the current results. 

Where possible, we have validated our results against quoted limits and generally find agreement 
within ~ 25%. Nonetheless there is the possibility of significant systematic errors since the detector 
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simulation we use does not fully capture the performance of the ATLAS and CMS detector, in 
particular in analyses that select for large missing energy and are therefore sensitive to the tails of 
distributions. For these reasons we emphasize that the results of this section should not be taken 
as strict experimental bounds but rather as approximate guidelines as to what regions of parameter 
space are accessible by various searches. 



We first consider constraints on the simplified model of section II C . We then discuss the realm of 
validity of the simplified model, considering further constraints on models with nonzero branching 
ratio of the NLSP to the Z and of the lightest squark directly to the LSP, as well as constraints 
on the decays of the heavier 3rd generation squarks. 



A. Constraints on the simplified model from SUSY searches 

CMS 3b + MET (a T ) 

The final states from our simplified model contain at least two 6-quarks, a moderate amount 
of missing energy, and the decay products of the two Higgses. A 125 GeV Higgs has a branching 
ratio to bottom quarks of 58% in the SM, so events can easily contain up to six 6-jets. Standard 
Model processes rarely produce more than two bottom quarks in an event, so searches in final 
states with three or more 6-jets can place strong limits on new physics. Both the ATLAS and CMS 
collaborations have performed searches in final states with three or more 6-tagged jets and missing 
energy [103] I104| . The CMS search in 3.9 fb _1 of 8 TeV data places the strongest individual 
limit; in the region Ht > 325 GeV, 11.7 events were expected and 6 were observed, for which 
we compute a 95% CL5 limit on the visible cross section of 1.3 fb. The results from our MC 
simulation for the visible cross section (total BSM cross section times branching ratio to final state 
times selection/identification efficiency) in this CMS search channel is shown in fig. [7J The bound 
is satisfied for NLSP masses < 200 GeV. The ATLAS 3b search places weaker bounds and is not 
shown. 

CMS SSDL + 2b 

Another type of analysis that can place relevant constraints are searches in final states with 6-jets 
and same-sign leptons, such as those performed by CMS in data at both 7 and 8 TeV (1051 1106] . 

This too is a rare final state in the Standard Model. In our models, same-sign leptons can be 
produced if both Higgses in the event produce a lepton in their decay. This rate is suppressed by 
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FIG. 7: CMS 3b + MET (cut): Visible cross section from our simplified model for the CMS search in 
final states with 3 6-jets and missing energy [104] . The thick dashed contour is the observed 95% limit, and 
dotted is the expected limit extrapolated to 20 fb -1 . 

small branching fractions, but unlike the 36-jet searches this analysis does not require large MET 
or Ht- The strongest constraint on our model comes from signal region 2 of the 8 TeV search |106| . 
which required two 6-tagged jets and two positively charged leptons as well as MET > 30 GeV 
and Ht > 80 GeV. Through an apparent downward fluctuation, zero events were observed in 
this channel though 11 were observed in the corresponding channel with two negatively charged 
leptons. A 95% CL upper limit of 2.8 expected events from non-SM processes was placed using 
the CLc, method, corresponding to a visible cross section of .71 fb at 8 TeV, which we interpret 

as a limit of 1.42 fb in the combined + + / channel (which would be populated evenly in our 

simplified model). The predicted cross sections from our simplified model are shown in fig. |8j For 
an NLSP mass < 160 GeV, the squark mass is bounded to be > 260 GeV. 

We validated our implementation of this search by computing efficiencies for the simplified 
model of pp — > bb, b — > txf,xf ~^ XiW^; the efficiencies we calculated corresponded to limits 20% 
stronger than reported in ref. [106J. 

CMS Multilepton 

Multilepton searches without 6-tags also put a strong constraint on our model, especially when 
identifying tau pairs, which are produced with branching fraction 6.4% in the decay of a 125 GeV 
Higgs. CMS has a performed search in three (3£) and four lepton (4£) channels in 4.98 fb -1 of 7 
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FIG. 8: CMS SSDL + 2b: Visible cross section from our simplified model for the cuts of signal region 
2 of the CMS search in final states with 2 6-jets and same-sign leptons at 8 TeV [106] . The thick dashed 
contour is the observed 95% limit, and dotted is the expected limit extrapolated to 20 fb -1 . 

TeV data |107| . Although our models have slightly higher rates in the 31 channels, much lower 
backgrounds make the 41 channels far more constraining. The most sensitive 41 channels are those 
vetoing Z resonances {'no-Z^), which are populated by our model when both Higgses decay either 
to WW* or r + r~. The search channels are further separated by the number of hadronic taus 
(0/1/2t) reconstructed. The branching ratios to the l/2r channels are several times larger than 
the Or channel while the background remains low. The dominant constraint over most of our 
parameter space comes from the 41, 2r, MET > 50 GeV, no-Z channels. Combining the low and 
high Ht channels, in 4.98 fb -1 2.2 ± 0.9 events were expected and 1 was observed, from which 
we compute a 95% CL5 limit on the visible cross section of 0.8 fb. Fig. [9] shows the visible cross 
section for our simplified model. The rate depends only weakly on the NLSP mass and the limit 
is satisfied for squark masses > 250 GeV. We note that the Or and lr regions with the same 
kinematic selections would also set a strong limit but for an observation of 4 events where ~ 1 
was expected. Taken together with the 2r channels, the observed rate would be consistent with a 
~ 250 GeV squark in our simplified model. We estimate the expected limit in 20 fb -1 of 8 TeV 
data (shown in fig. [6]) by applying the same search cuts to the 8 TeV signal while assuming the 
same background rate as in the 7 TeV search, scaled only by the luminosity. 

Refs. [108] consider the limits from the same multilepton search on SM hW production, and 
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FIG. 9: CMS Multilepton: Visible cross section from our simplified model for the combined 41, 2th, MET 
> 50 GeV, no-Z channels in the CMS 7 TeV multilepton search |107j . The thick dashed contour is the 
observed 95% limit. 

ref. [94J considers the limits set on hW production in BSM physics with a focus on the rare decay 
t — > ch. Compared to hW production, the multilepton limits on hh production are stronger due to 
the increased rates in 4£ channels and the importance of hadronic tau channels (not considered in 
refs. |94| I108j ). For comparison, our simulated efficiencies for the SM Higgs contribution to the Or 
channels were about 20% higher than reported in ref. [108] . likely due to a more detailed modeling 
of the CMS lepton efficiencies in their study. 

Limits in the 4.1 channels are statistics dominated and will likely improve significantly with the 
full 8 TeV run, especially with a properly combined analysis of the sensitive channels. If limits 
continue to improve, SM Higgs production will in fact be an important background in the most 
sensitive channels. The 3£ channels may also become sensitive if backgrounds can be reduced by 
6-tagging, a strategy also suggested in [9U [108J for identifying Higgs production in the tth and 
t —> ch channels. 

B. Constraints on NLSP branching ratio to Z 

An ATLAS search in 2.05 fb -1 for 6-jets and a leptonically decaying Z |109j places a constraint 
on the branching ratio of the NLSP to a Z in our model. This search targets pair production of 
stops in GMSB, with the stop decaying through the neutralino to a Z and gravitino, t — > b(xf — > 
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FIG. 10: 95% CL upper limit on the branching ratio of the NLSP to the Z in the simplified model from the 
ATLAS search in final states with b-jets and a leptonically decaying Z [109] . 

Xi) — > bZG. This is identical to our simplified model but with 100% branching fraction to a Z 
instead of to the Higgs. They quote a limit of > 310(330) GeV for an NLSP mass < (>) 190 
GeV. 

We determine efficiencies for this search by generating events in which one or both NLSPs 
decays to the Z instead of the Higgs. The visible cross section limits can be translated into a 95% 



constraint on the branching ratio BR(x!] — * Xi^Oj as shown in fig. 10 Branching ratios less than 



these values are easily achieved in the models described in sec. II A Nominally the contours of 



100% should agree with the quoted experimental exclusions for the GMSB stop model, and we see 
that our simulation is somewhat overestimating the strength of this search. 



C. Constraints on squark branching ratio directly to LSP 

If the squark can decay directly to the LSP instead of through the NLSP, then the bounds from 



searches with MET cuts and 6-tags become considerably tighter. Figure 11 show the 95% upper 
limit on the branching ratio to this mode from the CMS 36 + ot search, obtained by simulating 
events where one or both squarks decay directly to the LSP instead of through the NLSP. 

The limit on this branching ratio is 1-5% in the region of interest. For a weakly coupled LSP 
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FIG. 11: 95% CL upper limit on the branching ratio of the squark directly to the LSP from the CMS search 
in final states with 3 6-jets and missing energy in 3.9 fb _1 at 8 TeV [104J. The gray region is ruled out 
by this search even for 100% branching fraction to the NLSP (see fig. [7|. The dashed contour is the 95% 
exclusion for the specific model of a left handed sbottom with a wino NLSP and bino LSP. 



(photino/singlino/goldstino/gravitino) this constraint is easily satisfied. However, for a left handed 



sbottom with a wino NLSP and bino LSP, as described in sec. II B m- h < 350 GeV is ruled out. 



D. Constraints on heavier third generation squark decays 



While we have focused on the decays of the lightest squark, there are also constraints on the 
heavier squarks. When the lightest squark is bi, the constraints on directly constrain the 
parameter space of allowed sbottom masses (see eq. [TJ, while in the other cases the limits simply 
set a lower bound on the masses of other squarks in the spectra. 

The strongest constraints on the heavier squarks in the spectra occur when their decays give 
the final state tthh + MET, which can arise for example from the decay of a heavier stop to a bino 
NLSP, t —7- Bt with the B decaying to a Higgs. In particular the CMS SSDL + 2b search, discussed 
in sec. Ill A places a much stronger limit on this final state than on the bbhh + MET state due to 



leptonic decays of the W's produced from the top decay. We find that a squark decaying with 100% 
branching fraction to this mode with an on-shell top is excluded at 95% for masses < 400 GeV by 
the search in 3.9 fb -1 at 8 TeV; in 20 fb _1 the expected reach is ~ 450 GeV. 
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FIG. 12: Visible cross section for the decay thih — > W + W bb in the ATLAS single lepton direct stop search 
[110) . mi L and mg^ are related by eq. [TJ for which we chose tan/3 = 3. The thick dashed contour is the 
95% limit from 4.9 fb" 1 at 7 TeV. 



Direct stop searches have very targeted kinematic requirements and do not significantly con- 
strain the (tt/bb)hh + MET final states. However, in the case of ti decaying through a chargino 



co-NLSP to the final state W + W~bb + MET (as in figure 4(b)), the final state is similar enough 
to the ti + MET topology to be constrained. The strongest limit comes from the lowest MET 
signal region of the ATLAS 7 TeV search for stop pair production with a single leptonic W decay 
in 4.9 fb -1 [110J. When validating our method for the stop decay model in ref. |110j . we found 
the ratio of our efficiencies to those presented to be ~ 0.5 uniformly over the parameter space 
presented. We therefore scaled our efficiencies by a factor of 2 to compensate. For stops near 
350 GeV, the cross-section changes by a factor of two as the mass is changed by ~ 40 GeV, so we 
expect our bounds on squark mass to be within this margin of error. The constraints we obtain are 



shown in fig. 12 in terms of both the ti mass (top of frame) and the corresponding bi, mass (bot- 



tom of frame) for comparison with the direct constraints on bi production. The limit is roughly 
mr > 350 GeV for an NLSP heavier than 150 GeV. For heavier tr the tihh + MET final state 
also opens up, but we find that the dominant constraint is still from the direct stop search. 
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IV. DIPHOTON CHANNELS 



We now discuss the possibility of observing the simplified model of sec. II C in searches taking 
advantage of the Higgs diphoton resonance. We first describe our methods in the simplest case by 
considering the inclusive diphoton searches. We then discuss the potential reach of the existing 
CMS searches targeting SM associated production in 77 + t and 77 + MET channels and 
explore possible ways to increase their sensitivity to Higgs production in SUSY processes. Finally, 
we propose a search in 77 + 1/26 final states that could be much more sensitive to our model than 
other searches. 



A. Inclusive diphoton rates 

Both ATLAS and CMS have presented combined searches at 7 and 8 TeV for inclusive Higgs 
to diphoton final states, finding respectively best fit values of a/a^M = 1-9 ± 0.5 and ct/cxsm = 
1.56 ± 0.43 for a Higgs mass near 125 GeV [TT2l[TT3] . 

An important point is that this search does not veto on any extra activity in the event, so we 
expect our model to give a fractional enhancement to the rate of roughly ~ 2 x cr^/a^SMi where cr^ 
is the squark production cross section and the factor of 2 arises because two Higgses are produced 
in each SUSY event. To obtain a more accurate result reflecting the changes in kinematics and 
isolation between the SM and SUSY cascade production modes, we again generated simulated 



signals as described in sec. Ill and applied the kinematic cuts specified by each search. 



For searches in diphoton final states, the total efficiencies e depend on the efficiencies of the 
photon identification used by ATLAS and CMS. The PGS4 detector simulation applies only isolation 
cuts and a cut on the ratio of energy in the hadronic and EM calorimeters when identifying photons, 
while ATLAS and CMS use additional cuts on the shape and quality of the EM shower in the 
calorimeter. Our MC simulation therefore overestimates the photon identification efficiency. To 
account for this, we simulated pp —> Vh — > U77 in the SM and computed the efficiency for the 
ATLAS [112J and CMS |113j diphoton searches assuming the photon identification reported by 
PGS4. We define variables k as the ratios of these quantities (e]^~ ) to the actual pp — > Vh — > 
77 efficiencies reported by ATLAS and CMS including the tighter photon identification criteria 
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At both 7 and 8 TeV we find k atlas ~ ^cms ~ 0.7. We interpret k as a constant efficiency for 
both photons from a Higgs decay to pass the ATLAS or CMS photon ID criteria given that they 
pass the PGS4 photon ID criteria and the photon kinematic cuts, and normalize all efficiencies for 
diphoton searches computed using our MC simulation by this constant: 

,77 ~ ,77. MC fn\ 
BSM ~ t BSM x K \°) 

We chose to normalize to the Vh production mode because the Higgs has somewhat similar kine- 
matics to the BSM production modes we are considering. For comparison, applying the same 
procedure to the gg — > h production mode we obtain K gg -+h ~ 0.6. 

The rates from the search due to SM Higgs production (.Rsm) and SUSY Higgs production in 
our model (-Rbsm) are thus given by 

Rsm = o"sm x esM (4) 
A-Rbsm = 2crgg x £bsm = ^ a bb x e Bsi[ ( 5 ) 

where we assume the experimentally reported efficiencies and theoretically predicted cross section 
for a 125 GeV Higgs to determine the SM contribution. Note that we define e to include the 
branching ratio for the h — > 77 decay. 

Using the above method we can (approximately) determine the contribution of Higgs production 
in our simplified model to the total rate in the ATLAS and CMS inclusive h — )■ 77 searches. 
We define /i to be the ratio of the total signal rate including BSM production to the SM rate, 



^ = Rsm+aRbsm _ pjg_ shows fj, as a function of the squark mass and NLSP mass in the 8 TeV 
CMS inclusive diphoton search. The efficiencies are insensitive to the stop/sbottom and NLSP 
mass, and \x therefore depends dominantly on the squark production cross section. Because the 
squark cross section increases more rapidly with yfs then the Higgs cross section, at 7 TeV [i is 
about 10% smaller. The rates at ATLAS are similar to those at CMS. For a 250 GeV stop/sbottom, 
the total rate is ~ 150% of the SM rate. We will discuss the possibility of explaining current hints 
of excesses in the inclusive diphoton channel in section |V} 
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FIG. 13: (a): /j,, the ratio of total BSM + SM rate to the SM rate, for the CMS8 inclusive diphoton search, 
as a function of squark mass and NLSP mass in the simplified model. The shaded region indicates the 
existing 95% exclusion at ~ 5 fb _1 from SUSY searches and the diphoton + MET search, shown in detail 
m fig. § (b): The same quantity as a function of squark mass only, for a NLSP mass of 160 GeV. Shaded 
in green and blue are respectively the ATLAS and CMS la preferred regions. 

B. Diphoton plus lepton search 



CMS has published searches for diphoton resonances accompanied by a lepton or missing trans- 
verse energy (MET), targeting associated hW and hZ production with a focus on the fermiophobic 
Higgs model [111] . Both of these final states are promising search channels for observation of BSM 
Higgs production. 

In our simplified model, the signal in the 77+/ channel occurs when one Higgs decays to photons 
while the other produces leptons through decays to W^VF* , rr or ZZ* . Though the branching ratio 
to this state is small, the Standard Model background rate from non-Higgs processes and signal 
rate from hZ — > ^yj£ + £~ or hW — > 77^ production are both quite low. 

As before, we calculate the contribution to the rate in this search from our model, adjusting 



our MC efficiencies by kcms> and plot the results against squark mass. Fig. 14 shows the total rate 
both in units of cross section and as a ratio to the Standard Model rate. The search in 5.3 fb _1 of 8 
TeV data |lllj observed one event with about 1.6 expected background events within 125 ±2 GeV. 
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FIG. 14: The rate from our simplified model in the CMS diphoton plus lepton search |111) . both in events 
per fb and as a ratio [i to the rate from Standard Model Higgs production, as a function of squark mass for 
a NLSP mass of 160 GeV. The red shaded region is ruled out at 95% CL by the existing 5.3 fb -1 of data at 
8 TeV. Dashed is the expected limit in 20 fb -1 . 

With 0.5 SM Higgs events expected, this gives a 95% confidence level bound of 0.7 fb on the visible 
cross section from new physics, excluding squark masses < 240 GeV. Because the background is 
directly measured from the continuum diphoton spectrum, we scale the background uncertainty 
by yL to obtain an expected limit of 0.3 fb in 20 fb - . 

We validated our results by comparing our computed efficiencies for SM Wh and Zh events to 
the efficiencies determined from for a 120 GeV fermiophobic Higgs, finding agreement within 
10%. 

C. Diphoton plus missing energy search 

We now turn to the CMS search for diphoton + missing transverse energy (MET) targeting 
Zh — > vvjj |lllj . The Standard Model signal in this search is suppressed by branching ratios 
and the small Zh production cross section. However, this search can be quite sensitive to Higgs 
production in NLSP decays due to the guaranteed MET from the LSP, as long as the NLSP and 
Higgs are not too degenerate. 

The kinematic selections on the photons for this search are similar to the other diphoton chan- 
nels, except that only photons within the barrel are selected. Additionally, events passing the 
lepton tag are vetoed. Events with missing transverse energy greater than 70 GeV are selected. 
The background rate (fit from data) was 2.4 ± 0.8 events in the range 125 ± 2 GeV, and 3 events 
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were observed in 5.3 fb _1 . With a SM higgs expectation of 0.2 events, this gives a 95% limit on the 
visible cross section from new physics of 1.1 fb. The expected visible cross section for this search in 
the simplified model is shown as a function of squark and NLSP mass in fig. 15 In fig. [6| the 95% 



o 



exclusion are compared to the limits from SUSY searches, both for current data and extrapolated 
to 20 fb -1 . Since the background in this channel is fit from continuum diphoton data, we scaled the 
systematic uncertainty as yL to estimate the expected limits in 20 fb -1 . This search has better 
reach than the SUSY searches for NLSP masses of ~ 170 — 200 GeV. 

For NLSPs lighter than ~ 150 GeV, the signal in our model is suppressed significantly by a cut 
of MET > 70 GeV, which removes 80-90% of the signal events passing the other selections. The 
sensitivity of such a search can therefore likely be improved by including one or more additional 
channels with lower MET threshold, as was done in the non-resonant CMS diphoton + MET search 
[114J. As an example, the bin 40 GeV < MET < 70 GeV generally contains 2-3 times as many 
signal events as the MET > 70 GeV bin. If the background scales with the MET cut as in the 
non-resonant search [114] . it would be a factor of about 20 times larger in this bin, giving S/yB 
about half that of the higher MET channel. As the background increases rapidly with lower MET 
cuts, sub-dividing this region could further increase the sensitivity in lower MET channels. 

As before, we validated our results against the experimental efficiencies for a 120 GeV fermio- 
phobic Higgs. Our analysis yields a 125 GeV SM Higgs signal rate of R3^ MET = -04 fb, which 
agrees within 10% with the efficiency from the expected signal rate in [111] . 



D. Discovering SUSY in the (h — > 77) + 1/26 channel 

The Standard Model background rate for the h — I 77 search in 5.3 fb -1 at 8 TeV at CMS is 
433.8 events/GeV near the Higgs mass peak |113j . and consists dominantly of prompt diphoton 
(77) events and single prompt photon events with one or more jets (7J), with one jet faking a 
photon [115j . This background can be considerably reduced by requiring one or two 6-tags, while 
maintaining very good efficiency for the diphoton signal from our simplified model, which produces 
copious amounts of b quarks in the final state. We therefore consider a search with the same cuts 
as the CMS diphoton search, but with the additional requirement of one or two 6-tagged jets with 
Pt > 50 GeV. We consider a counting experiment in which the number of events with diphoton 
mass within 2 GeV of = 125 GeV is compared to the background expectation. The sensitivity 
of a similar search has been studied in refs. |116} 1117] in the context of measuring variations in the 
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FIG. 15: The visible cross section for our simplified model in the CMS diphoton plus MET search [121] as 
a function of squark mass and NLSP masses. We take the 95% limit from 5.3 fb -1 of data at 8 TeV to be 
1.1 fb. 

Higgs couplings in Higgs pair production events at very high luminosities with the 14 TeV LHC, 
and the sensitivity of a diphoton + jets search at 14 TeV without 6-tags has been considered in 
|31j for studying fermionic top-partner decays to the final states th + X. 

To estimate the continuum, non-Higgs background from 77 and 7J, we simulated yS = 8 TeV 
MLM matched samples of 77 + 0/l/2j and 7 + lj/2j/3j, where j is a light flavor or gluon, a 
charm quark, or b quark at 8 TeV. As before we simulate the response of the CMS detector with 
the default PGS CMS parameter card, but here we use a different modified 6-tagging function 
which gives a flat tag rate of 1% for light flavor/gluon jets, 15% for charm, and 60% for b quarks 
[118J. We separately normalize the 77 and 7J rates to the rate in |113j with relative proportions 
of 25% 7J and 75% 77 based on the Monte Carlo results in [115] . The results for the background 
rates are shown in table [TJ The first fo-tag reduces the non-Higgs background by a factor of 100 to 
fcont.,ife = 2.5 fb, with the remaining contributions from true 6-jets, light quark fakes, and charm 
fakes all roughly comparable. The second 6-tag reduces this background by another factor of 50 
to 0"cont.,2£> = 0.05 fb, with the remaining contributions dominantly from true 6-jets or mistagged 
charm. 

For the 26 selection, top pair production with hard radiated photons or electrons /jets faking 
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230 x 10 3 




130 200 


394 
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0.1 0.2 


0.2 
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16, p T > 50 GeV [fb] 
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0.6 


0.08 


2.5 


0.05 0.01 


0.01 


0.07 


1.4 


26, p T > 50 GeV [fb] 


0.025 0.007 


0.02 


0.05 


0.02 0.0004 0.0009 


0.02 


0.7 



TABLE I: The visible cross section for non-Higgs continuum and SM Higgs diphoton rates after various 
cuts in the 77 + 1/26 search, including all detector efficiencies, branching ratios, and fake rates. The 77 
and 7j categories include the non-resonant, non-electroweak contributions to the rate in events with 2 real 
photons or 1 real photon and 1 jet faking a photon respectively, tt is the contribution from top quark pair 
production with hard radiated photons or electrons /jets faking photons, tth is a Higgs-strahlung process in 
the SM, while bbh includes this and also QCD production of a 66 pair in gluon fusion Higgs production. For 
bbh we use the tree level cross section from MadGraph5. hZ is the associated production of a SM Higgs and 
Z —> bb. The cross section for tt is from [120| and for tth and hZ from |101j . The 77 and jj contributions 
are normalized to the observed background in [113] as described in the text. The signal rate in the simplified 
model for a 300 GeV squark and 150 GeV NLSP is also shown. 



photons is also an important continuum background. To estimate this, we take a jet-photon fake 
rate of 0.1% |119] and an electron-photon fake rate of 2% |114j . The estimated rate is cr t j 2b = 
0.02 fb, with fake electrons and radiated photons giving the dominant contributions. SM tth 
also gives a resonant contribution to the background with a rate comparable to the continuum 
background. Zh and bbh production give a negligible rate, as does non-resonant electroweak 
production with prompt or fake photons. The 26 channel is therefore essentially backgroundless 
at 5 fb -1 and in 20 fb -1 has ~ 1 expected continuum event falling in the Higgs mass window and 
< 1 expected SM tth event. Importantly, because this is a search for a known resonance, the 
background rate can be measured from data as in the usual Higgs diphoton search, so an actual 
search would be limited only by statistical uncertainties in the background fit. 

The visible cross section for the lb and 2b selections are shown for our simplified model in fig. 
16 The efficiency for signal events passing the inclusive diphoton selections to also pass the lb 
selection is 50 — 70%, while for the 2b selection it is 20 — 40%. For the lb channel we estimate the 
systematic uncertainty in a background fit to be 20% and 10% at 5 fb -1 and 20 fb -1 respectively 
based on counting statistics in the neighboring diphoton invariant mass bins. With a background 
rate of 0bg,i& = 2.5 fb, we find expected 95% CL S limits on the visible cross of 1.9 fb for 5 fb -1 
and 0.9 fb -1 for 20 fb -1 . The 16 channel has sensitivity significantly exceeding the other diphoton 
and SUSY searches we have discussed, as shown in fig. [6j 
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FIG. 16: The visible cross section for our simplified model in a search for h —¥ 77 in association with 
6-tagged jets, (a) Requiring at least one 6-tagged jet with px > 50 GeV. (b) Requiring at least two 6-tagged 
jets with pt > 50 GeV . 

The near backgroundless 26 channel is even more sensitive to our model and can even lead 
to a discovery at the 8 TeV LHC. Since there are large uncertainties in our background estimate 
due to the fake rates for photons and 6-tags, we estimate the reach of the 26 channel taking the 



background rate to be 0.1±0.05 fb. Fig. 17 shows the expected signal significance in 5 and 20 fb -1 , 
along with the 95% CL excluded region from the existing SUSY and diphoton+MET searches and 
their extrapolated limits in 20 fb -1 . For direct comparison with the other searches, the expected 
95% exclusions from this channel are also shown in fig. |6| 

In 5 fb -1 the 77 + 26 search already has significantly better reach than the other searches we 
have considered, and in 20 fb -1 it gives an expected significance of > 3a for squarks up to 350 GeV, 
well beyond even the exclusion reach for the other searches, and a potential 5a discovery for squarks 
up to 300 GeV. Another important advantage of this channel over other SUSY searches that can 
probe our models is that by exploiting the Higgs diphoton resonance, this search clearly identifies 
the the involvement of an on-shell Higgs boson in the new physics. In contrast, a signal observed 
in the multi-6 or leptonic searches would not immediately suggest excess Higgs production. 
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FIG. 17: The expected signal significance for a proposed search in 77 + 26 in 8 TeV with (a) 5 fb" 1 and (b) 
20 fb" 1 . The shaded region indicates (a) the existing 95% exclusions at ~ 5 fb 1 and (b) the extrapolated 
limits at 20 fb" 1 from SUSY searches and the diphoton + MET search, as described in fig. [6j 

V. DIPHOTON EXCESSES AND OTHER LHC HIGGS SEARCHES 



Ref. [9] presents a combined fit to the ATLAS and CMS 7 and 8 TeV signals in the various SM 
Higgs searches. The only > la deviations from the SM predictions are excesses in the inclusive 
and VBF tagged 77 channels and a low rate in the rr channel, driven primarily by a low rate in 
the VBF tagged h -> rr channel at CMS [2]. 



As shown in fig. 13 for a squark mass of 250-290 GeV, an enhancement of 20-50% in the 
diphoton channel from Higgs production in SUSY cascades is possible while remaining consistent 
with other searches. This search is sensitive to our model because it makes use of the clean 77 
resonance and does not place cuts on the non-Higgs activity in the event. This is also true for the 
h —7- ZZ* — > At resonance searches, so we would expect any enhancement in the inclusive diphoton 
searches by this mechanism to give a comparable enhancement to h — > ZZ* — > A.I. 

On the other hand, searches in other channels often require additional cuts on the event to 
reduce SM backgrounds and consider more complicated kinematic variables, resulting in lowered 
signal strength from the non-standard Higgs production modes we consider. For example, in order 
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to reduce backgrounds from top production, the ATLAS and CMS searches in the h — > WW* — > 
lulu \\.22\ 1123] and in most channels for h — )■ tt (1241 11251 veto events with any 6-tagged jets. 
Since events from our models typically contain several 6-jets, they give a negligible contribution to 
the rate in these searches. Furthermore, the kinematic variables considered in these analyses (such 
as the transverse mass of leptons and the reconstructed tau-tau invariant mass) are distorted by 
the different kinematics and additional missing energy in events from squark production, so those 
events that are selected will not necessarily increase the signal strength for the 125 GeV SM Higgs 
hypothesis. Similar issues limit the sensitivity of searches for Vh production with h — > bb and the 
vector boson V decaying to leptons and/or neutrinos |126| il27j. Although our model produces 
these same final states, the 6-jets selected by the analyses will not necessarily be the Higgs decay 
products, so a mass peak will not be visible. Both analyses place many kinematic cuts on events, 
such as requiring the bb momentum to be nearly back-to-back with the reconstructed vector boson 
momentum when projected into the transverse plane, greatly reducing the efficiency in our case. 

In the h — >■ 77 analyses both CMS and ATLAS define dijet-tagged event classes, selecting jets 
widely separated in pseudorapidity r] to obtain sensitivity to VBF Higgs production in the SM. The 
additional jets present in BSM events from our simplified model rarely satisfy this Ar/ requirement, 
and the enhancement relative to the SM rate is at most a few percent in the parameter range we 
consider. 

To summarize, if the simplified model of excess Higgs production we have considered gives an 
observable enhancement to inclusive h — >■ 77 rate, then this is consistent with the lack of observed 
excesses in h — > WW* , h — > tt, and Vh —> Vbb, but is in slight tension with an observed rate for 
h — > Ai which is consistent with the SM and does not explain the excess in the dijet-tagged h — > 77 
channel. Therefore, while the enhancement to the inclusive diphoton channel is encouraging, the 
simplified model alone likely does not give a significantly better fit to current data than does the 
SM prediction. 

Thus far we have always assumed that the effective couplings of the Higgs are the same as those 
predicted by the Standard Model. However, the light third generation squarks of our models can 
modify the Higgs couplings to gluons and photons through loop effects, particularly if there are light 
stops (discussed recently in refs. [8, 128-130J). Large stop mixing (e.g. from ^4-terms) is required 
to give the correct sign to enhance the diphoton branching ratios, and extremely large mixing is 
required to actually enhance the total diphoton rate from gluon fusion Higgs production since the 
gg — >■ h rate is at first suppressed as the diphoton branching ratio is enhanced [8] . By contributing 
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another production mode which is not sensitive to the gg — )• h coupling, Higgs production in cascade 
decays allows an overall enhancement to the diphoton rate at smaller mixings. The modified h — > 77 
branching ratio then also enhances the dijet-tagged diphoton rate to reproduce the observed excess, 
while the reduced gg — )• h rate prevents an excess in the h — > ZZ* — > 4/ channel. 

We do not attempt a detailed fit to the data, but rather note that in the range of squark 
masses where these models can give significant enhancements to inclusive SM Higgs searches, it 



is likely that excesses in the 77 + MET search and SUSY searches discussed in sec. IV C and III 



will be observed in 20 fb 1 of 8 TeV data and that a discovery could be made in a search for 



77 + 1/26 as described in sec. IV D As a general point, excess Higgs production in BSM cascades, 
possibly with less dramatic exclusive signatures than in our models, is an interesting degree of 
freedom to consider when fitting rates in Higgs searches, giving qualitatively different effects from 
modifications to Higgs couplings. 

We also note that to improve sensitivity, both ATLAS and CMS divide the inclusive diphoton 
search into several event classes based on the kinematics of the photon and the quality of its 



reconstruction. Our results in sec. IV A based on the inclusive diphoton count therefore only give 
an approximation of expected best fit rates, which in refs. [1121 1113] are actually determined with 
a fit taking into account the different signal and background expectations in the individual event 
classes. Of particular note, both ATLAS and CMS have event classes capturing Higgs production 
with higher transverse momentum, which has a smaller background rate. For example, ATLAS 
separates diphoton candidates into high and low pxt event classes, where pxt is the component of 
the diphoton momentum transverse to the thrust axis [112j . The simplified model Higgs production 
leads to a substantially larger signal enhancement in the high prt category compared to SM gg — > h 
events. For a given enhancement to the inclusive rate relative to the SM, the expected enhancement 
in the high pxt categories relative to the SM rate in these channels is generally 4 — 5 times greater. 
In the combined 7 + 8 TeV results these channel show roughly equal enhancements compared to 
the SM rate, but with large enough individual uncertainty to accommodate a substantially higher 
rate in the high pj-t channels. While other exclusive diphoton channels are more sensitive to our 
specific model, with more statistics this could be a useful constraint on more general possibilities 
for enhancing Higgs production through BSM cascades. 
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VI. DIRECT PRODUCTION OF NEUTRALINOS AND CHARGINOS 

If the NLSP in our models is a sufficiently light wino or higgsino, then the electroweak direct 
pair production rates can be large compared to SM associated Higgs production rates. This can 
mildly constrain the neutralino/chargino spectra possible for the colored production case, and can 
also enhance the signal in diphoton searches. As discussed in sec. |II A[ we consider the possibilities 
of wino co-NLSPs with small splittings such that the charginos decay promptly to the LSP, and a 
mixed wino-bino or higgsino NLSP, with the nearby neutralinos and charginos decaying promptly 
to the NLSP. 

For wino-like co-NLSPs, the dominant production modes are X2X1 an d X\Xi (here X2 denotes 
the neutral NLSP and x? the LSP). The winos decay as Xi ~^ W^xi an< ^ X2 hxi, so the final 
states of interest are hW + MET and W + W~ + MET. Refs. [133 EM] discuss the possibility that 
X\Xi production can appear in h — > W + W~ search channels, but for X\ decaying dominantly to 
there is no significant constraint on direct production in the parameter space of interest. The 
discovery potential for heavier charginos and neutralinos in the hW mode at the 8 and 14 TeV 
LHC has been studied in refs. |33| 134] . However for a lighter NLSP with small splitting with the 
Higgs, the MET signature is suppressed and the resulting signals become much more similar to 
SM associated Higgs production. Because the production cross section can be many times larger 
than the SM rate and there is no additional activity in the event (unlike in the colored production 
case), all of the existing searches for SM-like hW associated production will directly constrain the 
X2X1 production cross section, which is shown relative to the SM hW cross section at 8 TeV in 
fig. [l8j The most constraining direct limit on hW is an ATLAS search is for hW — > 3l3v at 7 TeV 
with a limit of crwh/&Wh,SM ^ 6 [135] . requiring > 135 GeV; a similar bound in the same 
channel is also presented in ref. [136J in a search for direct production of gauginos with leptophilic 
decays. The limit from the CMS 77 + I search [111] discussed earlier is o"wh/cwh,sm ^ 7. The 
strongest combined constraint on Zh/Wh production currently comes from a CMS search for 
h — > bb with leptonically decaying VF's and Z's at 7 and 8 TeV, which gives a 95% upper limit of 
o"vh/cvh,sm < 2.11 for a 125 GeV Higgs [137] . but individual bounds on hW and hZ production 
were not presented. Ref. [35] has suggested using jet substructure methods to further improve 
sensitivity for direct production of gauginos decaying in the bbW channel. 

For a mixed wino-bino NLSP or a higgsino-like NLSP, the charginos can promptly decay to the 
NLSP which then decays to a weakly coupled LSP, so the dominant signatures are in (hh/hZ/ZZ) + 
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W + W° -» hW + ~+MET 
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FIG. 18: The cross section for wino co-NLSP direct production relative to the SM cross section for hW 
associated production at 8 TeV. Dashed is the 95% CL limit on the BSM contribution in the ATLAS 
hW — > 3l3is channel, and dotted is the limit from the CMS 77 + I channel. 
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FIG. 19: Rates in diphoton channels for direct production. Solid is a mixed wino-bino NLSP, dotted is 
a wino co-NLSP, dashed is a Higgsino NLSP at tan/3 = 1.5, and dashed dotted is a Higgsino NLSP at 
tan/3 = 1.1. (a) Inclusive CMS8 diphoton rate relative to the SM Higgs rate, (b) Visible cross section in 
CMS8 77 + MET > 70 GeV search. Dashed horizontal is the expected 95% limit in 20 fb" 1 . (c) Visible 



cross section in 77 + 16 search described in sec. IV D Dashed horizontal is the expected 95% limit in 20 fb 
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MET. The Tevatron and LHC signatures for direct production of higgsinos decaying to gravitinos 
have been discussed in [38H40| . with emphasis on the situations where the branching ratio to Z 
dominates or is comparable to the Higgs branching ratio. We focus instead on the case where the 
hh mode is dominant, which in the higgsino case at low masses requires tan/3 ~ 1. These modes 
are even less constrained than the wino co-NLSP because of their smaller branching fraction to 
leptonic final states. Refs. |38[I39| found that a dedicated Tevatron search in hh — > 46+ MET could 
be sensitive to a higgsino NLSP with masses ~ 150 — 200 GeV, but to date the most constraining 
search has been a CDF search |138j for gluinos decaying to sbottoms requiring 26-tags and MET > 
70 GeV, which does not constrain either case even assuming the same efficiencies as the targeted 
g —> bb channel. Ref. [ID] briefly discussed the possibility of a 46 + MET search at the LHC, but 
again we find that existing searches (CMS Qy and ATLAS 36 + MET) do not set a limit on direct 
production. 

The signals in the inclusive diphoton, diphoton+MET, and diphoton+16 channels are shown in 



fig. 19 for the two wino NLSP possibilities assuming 100% branching fractions for the NLSP to the 
Higgs, and for a higgsino NLSP decaying to a gravitino LSP at tan/3 = 1.1 and tan/3 = 1.5. Only 
the mixed wino-bino NLSP with prompt decays of the charginos to the NLSP has a potentially 
interesting signal in these channels, due to its large production cross section compared to a higgsino 
NLSP and the potential to produce two Higgses in each event, unlike the wino co-NLSP scenario. 
For masses lighter than 150 GeV it can have a signal of 20-30% the SM rate in the inclusive 
77 channel and can be potentially observed in the 77 + 16 channel in 20 fb . The rate in the 
77 + 26 channel is too low to afford sensitivity. The 77 + MET search is just barely sensitive 
enough to excludes masses around 200 GeV in 20 fb -1 , but for lower masses the MET signal is too 
suppressed. Again, the rate in a lower MET bin can be several times higher and may increase the 
sensitivity for a lighter wino. Although the rates for the other NLSP scenarios are too small to be 
observable from just direct production, in the context of the simplified squark production model 
the additional contribution from direct production can slightly enhance the overall sensitivity of 
the diphoton channels. 



VII. CONCLUSIONS 

We have shown that in a variety of models the decays of a neutralino NLSP can give copious 
extra Higgs boson production while remaining consistent with constraints from BSM searches. This 
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leads to potentially observable signatures in Higgs resonance searches. Well- motivated realizations 
of this scenario include wino- and bino-like NLSPs with a singlino/goldstino/string photino LSP, 
or a higgsino-like NLSP with a gravitino LSP. We have focused on pair production of light third 
generation squarks, motivated by naturalness of the electroweak scale. A wide class of realistic 
models can be well described by the simplified model of fig. [5j leading to final states with two 
bottom quarks, two Higgses, and missing energy. Existing searches constrain these models to have 
squark masses > 260 GeV and NLSP masses < 200 GeV, as summarized in fig. [6} For light squarks 
we argue that the most effective search strategy to probe these models is to target the h — > 77 decay 
mode, which would give a clear signal of a Higgs resonance. A search in 77 + 2b final states would 
have high signal efficiency and very low SM background, and could give high signal significance in 



8 TeV LHC data for squarks up to ~ 350 GeV (figure 17). 

Although a 77 + 2b search is expected to be the best probe for the models we consider, even 
existing analyses can have some sensitivity with further 8 TeV data. Of particular interest is the 
possibility that the BSM production rate is high enough to account for the high signal strength 
reported by ATLAS and CMS in the inclusive diphoton channels. This would require squarks near 



the limits placed by SUSY searches (see fig. 13). Some excess are then expected for these searches 
in 8 TeV data, though this would be difficult to link to Higgs physics. More informative would be 
the searches for associated Higgs production, namely the 77 + £ and 77 + MET analyses, which 
would all have rates well in excess of the Standard Model prediction. As discussed in sec. |VJ little 
or no enhancement is expected in most searches in other Higgs decay modes such as WW* or 
rr. Therefore, observing excesses in only the diphoton channels and particularly the associated 
production channels would suggest our class of models. In the parameter space where observable 
excesses can occur in these channels, a search in 77 + 2b can likely give a 5a discovery in 20 fb . 

Although we have focused on a particular SUSY simplified model in this work, similar phe- 
nomenology can arise in many other models. Final states very similar to those in our simplified 
model can arise for example from decays of new vector-like quarks to a Higgs and third generation 
quarks, especially in the decays of bottom-like quarks, e.g. b' —> hb. More generally, much of what 
we have said is relevant for pair production of any new particles that decay dominantly to Higgses. 
The diphoton decay mode still offers an excellent opportunity for identifying the Higgs resonance, 
and a 77 + lb search can still have good efficiency because of the large Higgs branching ratio to bb. 

Finally, while this work has explored one well motivated scenario in which non-standard Higgs 
channels can become sensitive, we emphasize the importance of searching for Higgs production in 
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as many channels as possible. Abandoning all theoretical prejudice, little is currently known about 
the production and decay of the newly observed resonance at 125 GeV, and there remains open 
the possibility of observing significant departures from Standard Model expectations in the very 
near future. 

Acknowledgments 

We would like to thank Masha Baryakhtar, Timothy Cohen, Nathaniel Craig, Savas Dimopou- 
los, Peter Graham, Jo-Anne Hewett, Caitlin Malone, John March-Russell, Tom Rizzo, Surjeet 
Rajendran, Ken van Tilburg and Giovanni Villadoro for useful discussions. We especially thank 
Timothy Cohen and Peter Graham for giving comments on a draft of this work before publication. 
P.S. is supported by a Stanford Graduate Fellowship, and K.H. is supported by a NSF Graduate 
Research Fellowship under Grant No. DGE-0645962. This work was partially supported by ERC 
grant BSMOXFORD no. 228169. 



[1] G. Aad et al. [ATLAS Collaboration], arXiv:1207.7214 [hep-ex]. 

[2] [CMS Collaboration] Evidence for a new state decaying into two photons in the search for the standard 

model Higgs boson in pp collisions. CMS-PAS-HIG-12-015 [2012] 

[3] A. Azatov, R. Contino and J. Galloway, JHEP 1204, 127 (2012) [arXiv: 1202.3415 [hcp-ph]]. 

[4] R. Bcnbrik, M. G. Bock, S. Heinemeyer, O. Stal, G. Weiglein and L. Zeune, arXiv:1207.1096 [hcp-ph]. 

[5] D. Carmi, A. Falkowski, E. Kuflik, T. Volansky and J. Zupan, arXiv:1207.1718 [hep-ph]. 

[6] J. Ellis and T. You, JHEP 1206, 140 (2012) [arXiv:1204.0464 [hcp-ph]]. 

[7] J. R. Espinosa, C. Grojean, M. Muhlleitner and M. Trott, arXiv:1207.1717 [hep-ph]. 

[8] P. P. Giardino, K. Kannikc, M. Raidal and A. Strumia, JHEP 1206, 117 (2012) [arXiv: 1203.4254 

[hep-ph]]. 

[9] P. P. Giardino, K. Kannike, M. Raidal and A. Strumia, arXiv:1207.1347 [hep-ph]. 

[10] M. Klute, R. Lafayc, T. Plehn, M. Rauch and D. Zcrwas, arXiv:1205.2699 [hep-ph]. 

[11] I. Low, J. Lykkcn and G. Shaughnessy, arXiv: 1207. 1093 [hcp-ph]. 

[12] F. del Aguila, G. L. Kane and M. Quiros, Phys. Rev. Lett. 63, 942 (1989). 

[13] F. del Aguila, L. Ametller, G. L. Kane and J. Vidal, Nucl. Phys. B 334, 1 (1990). 

[14] M. Shcr, Phys. Rev. D 61, 057303 (2000) [hcp-ph/9908238]. 

[15] J. A. Aguilar-Saavcdra, JHEP 0612, 033 (2006) [hcp-ph/0603200]. 

[16] C. Kilic, K. Kopp and T. Okui, Phys. Rev. D 83, 015006 (2011) [arXiv: 1008.2763 [hcp-ph]]. 



38 



[17] A. Carmona, M. Chala and J. Santiago, JHEP 1207, 049 (2012) [arXiv: 1205.2378 [hep-ph]]. 

[18] G. D. Kribs, A. Martin, T. S. Roy and M. Spannowsky, Phys. Rev. D 81, 111501 (2010) 

[arXiv:0912.4731 [hep-ph]]. 
[19] P. Bandyopadhyay, E. J. Chun and J. -C. Park, JHEP 1106, 129 (2011) [arXiv:1105.1652 [hep-ph]]. 
[20] M. Son, C. Spethmann and B. Tweedie, arXiv: 1204.0525 [hep-ph]. 

[21] P. Bandyopadhyay,Bandyopadhyay:2010tv JHEP 1108, 016 (2011) [arXiv: 1008.3339 [hep-ph]]. 
[22] A. Datta, A. Djouadi, M. Guchait and Y. Mambrini, Phys. Rev. D 65, 015007 (2002) [hcp-ph/0107271]. 
[23] A. Datta, A. Djouadi, M. Guchait and F. Moortgat, Nucl. Phys. B 681, 31 (2004) [hcp-ph/0303095]. 
[24] F. dc Campos, O. J. P. Eboli, M. B. Magro, D. Restrcpo and J. W. F. Valle, Phys. Rev. D 80, 015002 

(2009) [arXiv:0809.1637 [hep-ph]]. 
[25] A. C. Fowler and G. Weiglein, JHEP 1001, 108 (2010) [arXiv:0909.5165 [hep-ph]]. 
[26] S. Gori, P. Schwallcr and C. E. M. Wagner, Phys. Rev. D 83, 115022 (2011) [arXiv:1103.4138 [hep-ph]]. 
[27] G. D. Kribs, A. Martin, T. S. Roy and M. Spannowsky, Phys. Rev. D 82, 095012 (2010) 

[arXiv:1006.1656 [hep-ph]]. 
[28] O. Stal and G. Weiglein, JHEP 1201, 071 (2012) [arXiv:1108.0595 [hep-ph]]. 

[29] G. Azuelos, K. Benslama, D. Costanzo, G. Couture, J. E. Garcia, I. Hinchliffc, N. Kanaya and 

M. Lechowski et al, Eur. Phys. J. C 39S2, 13 (2005) [hcp-ph/0402037]. 
[30] G. D. Kribs, A. Martin and T. S. Roy, Phys. Rev. D 84, 095024 (2011) [arXiv:1012.2866 [hep-ph]]. 
[31] A. Azatov, O. Bondu, A. Falkowski, M. Felcini, S. Gascon-Shotkin, D. K. Ghosh, G. Moreau and 

S. Sekmcn, arXiv:1204.0455 [hep-ph]. 
[32] N. Vignaroli, JHEP 1207, 158 (2012) [arXiv:1204.0468 [hep-ph]]. 

[33] H. Baer, V. Barger, A. Lessa, W. Sreethawong and X. Tata, Phys. Rev. D 85, 055022 (2012) 

[arXiv:1201.2949 [hep-ph]]. 
[34] D. Ghosh, M. Guchait and D. Sengupta, arXiv:1202.4937 [hep-ph]. 
[35] P. Byakti and D. Ghosh, arXiv:1204.0415 [hep-ph]. 

[36] M. Asano, H. D. Kim, R. Kitano and Y. Shimizu, JHEP 1012, 019 (2010) [arXiv:1010.0692 [hep-ph]]. 
[37] Y. Kats, P. Meade, M. Reece and D. Shih, JHEP 1202, 115 (2012) [arXiv:1110.6444 [hep-ph]]. 
[38] K. T. Matchev and S. D. Thomas, Phys. Rev. D 62, 077702 (2000) [hep-ph/9908482]. 
[39] P. Meade, M. Reece and D. Shih, JHEP 1005, 105 (2010) [arXiv:0911.4130 [hep-ph]]. 
[40] J. T. Rudcrman and D. Shih, arXiv: 1103.6083 [hep-ph]. 

[41] J. Thaler and Z. Thomas, JHEP 1107, 060 (2011) [arXiv:1103.1631 [hep-ph]]. 

[42] D. Das, U. Ellwanger and A. M. Tcixcira, JHEP 1204, 067 (2012) [arXiv: 1202.5244 [hep-ph]]. 

[43] M. Baryakhtar, N. Craig and K. Van Tilburg, arXiv:1206.0751 [hep-ph]. 

[44] J. A. Aguilar-Saavedra, PoS TOP 2006, 003 (2006) [hcp-ph/0603199]. 

[45] J. A. Aguilar-Saavedra, JHEP 0911, 030 (2009) [arXiv:0907.3155 [hep-ph]]. 

[46] K. Harigaya, S. Matsumoto, M. M. Nojiri and K. Tobioka, arXiv:1204.2317 [hep-ph]. 

[47] A. Girdhar and B. Mukhopadhyaya, arXiv:1204.2885 [hep-ph]. 



39 



[48] N. Vignaroli, arXiv: 1207.0830 [hep-ph]. 

[49] H. Baer, M. Bisset, X. Tata and J. Woodside, Phys. Rev. D 46, 303 (1992). 

[50] P. Bandyopadhyay, A. Datta and B. Mukhopadhyaya, Phys. Lett. B 670, 5 (2008) [arXiv:0806.2367 
[hep-ph]]. 

[51] P. Bandyopadhyay, JHEP 0907, 102 (2009) [arXiv:0811.2537 [hep-ph]]. 
[52] B. Bhattachcrjec and A. Datta, JHEP 1203, 006 (2012) [arXiv: 1104.4529 [hep-ph]]. 
[53] I. Hinchliffc, F. E. Paige, M. D. Shapiro, J. Soderqvist and W. Yao, Phys. Rev. D 55, 5520 (1997) 
[hep-ph/9610544]. 

[54] K. Huitu, R. Kinnunen, J. Laamanen, S. Lehti, S. Roy and T. Salminen, Eur. Phys. J. C 58, 591 

(2008) [arXiv:0808.3094 [hep-ph]]. 
[55] [ATLAS Collaboration], "Search for squarks and gluinos with the ATLAS detector using final states 

with jets and missing transverse momentum and 4.7 fb _1 of ^/s = 7 TcV proton-proton collision 

data," ATLAS-CONF-2012-033 (2012). 
[56] [CMS Collaboration], "Interpretation of Searches for Supersymmetry" CMS-PAS-SUS-11-016 (2011). 
[57] S. Dimopoulos and G. F. Giudice, Phys. Lett. B 357, 573 (1995) [hep-ph/9507282]. 
[58] A. Pomarol and D. Tommasini, Nucl. Phys. B 466, 3 (1996) [hep-ph/9507462]. 
[59] A. G. Cohen, D. B. Kaplan and A. E. Nelson, Phys. Lett. B 388, 588 (1996) [hep-ph/9607394]. 
[60] K. Agashe and M. Graesser, Phys. Rev. D 59, 015007 (1999) [hcp-ph/9801446]. 

[61] D. E. Kaplan, F. Lepcintre, A. Masiero, A. E. Nelson and A. Riotto, Phys. Rev. D 60, 055003 (1999) 
[hep-ph/9806430]. 

[62] O. Aharony, L. Berdichevsky, M. Bcrkooz, Y. Hochberg and D. Roblcs-Llana, Phys. Rev. D 81, 085006 

(2010) [arXiv: 1001.0637 [hep-ph]]. 
[63] R. Barbieri, E. Bertuzzo, M. Farina, P. Lodonc and D. Zhuridov, JHEP 1012, 070 (2010) [Erratum- 

ibid. 1102, 044 (2011)] [arXiv: 101 1.0730 [hep-ph]]. 
[64] N. Craig, D. Green and A. Katz, JHEP 1107, 045 (2011) [arXiv: 1103.3708 [hep-ph]]. 
[65] G. Larsen, Y. Nomura and H. L. L. Roberts, JHEP 1206, 032 (2012) [arXiv:1202.6339 [hep-ph]]. 
[66] N. Craig, S. Dimopoulos and T. Ghcrghctta, JHEP 1204, 116 (2012) [arXiv:1203.0572 [hep-ph]]. 
[67] N. Craig, M. McCullough and J. Thaler, JHEP 1206, 046 (2012) [arXiv:1203.1622 [hep-ph]]. 
[68] H. Baer, V. Bargcr, P. Huang and X. Tata, JHEP 1205, 109 (2012) [arXiv:1203.5539 [hep-ph]]. 
[69] T. Cohen, A. Hook and G. Torroba, arXiv:1204.1337 [hep-ph]. 
[70] L. Randall and M. Rcecc, arXiv:1206.6540 [hep-ph]. 

[71] R. Dcrmisck, Mod. Phys. Lett. A 24, 1631 (2009) [arXiv:0907.0297 [hep-ph]]. 

[72] P. W. Graham, A. Ismail, S. Rajendran and P. Saraswat, Phys. Rev. D 81, 055016 (2010) 

[arXiv:0910.3020 [hep-ph]]. 
[73] S. P. Martin, Phys. Rev. D 81, 035004 (2010) [arXiv:0910.2732 [hep-ph]]. 

[74] K. Choi, K. S. Jcong, T. Kobayashi and K. i. Okumura, Phys. Lett. B 633, 355 (2006) [arXiv:hep- 
ph/0508029]. 



40 



[75] R. Kitano and Y. Nomura, Phys. Lett. B 631, 58 (2005) [arXiv:hep-ph/0509039]. 
[76] Z. Chacko, Y. Nomura and D. Tuckcr-Smith, Nucl. Phys. B 725, 207 (2005) [arXiv:hep-ph/0504095]. 
[77] J. R. Ellis, J. F. Gunion, H. E. Habcr, L. Roszkowski and F. Zwirner, Phys. Rev. D 39, 844 (1989). 
[78] J. R. Espinosa and M. Quiros, Phys. Rev. Lett. 81, 516 (1998) [arXiv:hep-ph/9804235]. 
[79] P. Batra, A. Delgado, D. E. Kaplan and T. M. P. Tait, JHEP 0402, 043 (2004) [arXiv:hcp-ph/0309149]. 
[80] A. Maloney, A. Pierce and J. G. Wacker, JHEP 0606, 034 (2006) [arXiv:hcp-ph/0409127]. 
[81] J. A. Casas, J. R. Espinosa and I. Hidalgo, JHEP 0401, 008 (2004) [arXiv:hcp-ph/0310137]. 
[82] A. Brignole, J. A. Casas, J. R. Espinosa and I. Navarro, Nucl. Phys. B 666, 105 (2003) [arXiv:hep- 
ph/0301121]. 

[83] R. Harnik, G. D. Kribs, D. T. Larson and H. Murayama, Phys. Rev. D 70, 015002 (2004) [arXiv:hep- 
ph/0311349]. 

[84] S. Chang, C. Kilic and R. Mahbubani, Phys. Rev. D 71, 015003 (2005) [arXiv:hep-ph/0405267]. 
[85] A. Delgado and T. M. P. Tait, JHEP 0507, 023 (2005) [arXiv:hcp-ph/0504224]. 
[86] A. Birkedal, Z. Chacko and Y. Nomura, Phys. Rev. D 71, 015006 (2005) [arXiv:hep-ph/0408329]. 
[87] K. S. Babu, I. Gogoladze and C. Kolda, arXiv:hcp-ph/0410085. 

[88] K. Choi, K. S. Jeong, T. Kobayashi and K. i. Okumura, Phys. Rev. D 75, 095012 (2007) [arXiv:hep- 
ph/0612258]. 

[89] C. Brust, A. Katz, S. Lawrence and R. Sundrum, JHEP 1203, 103 (2012) [arXiv: 11 10.6670 [hcp-ph]]. 

[90] M. Papucci, J. T. Ruderman and A. Weiler, arXiv:1110.6926 [hep-ph]. 

[91] H. M. Lee, V. Sanz and M. Trott, JHEP 1205, 139 (2012) [arXiv:1204.0802 [hep-ph]]. 

[92] A. Arvanitaki, N. Craig, S. Dimopoulos, S. Dubovsky and J. March-Russell, Phys. Rev. D 81, 075018 

(2010) [arXiv: 0909. 5440 [hep-ph]]. 
[93] C. Brust, A. Katz and R. Sundrum, arXiv:1206.2353 [hep-ph]. 

[94] N. Craig, J. A. Evans, R. Gray, M. Park, S. Somalwar, S. Thomas and M. Walker, arXiv: 1207.6794 
[hep-ph]. 

[95] G. D. Kribs and A. Martin, arXiv: 1207.4496 [hep-ph]. 

[96] J. Alwall, M. Herquct, F. Maltoni, O. Mattelaer and T. Stelzer, JHEP 1106, 128 (2011) 
[arXiv:1106.0522 [hep-ph]]. 

[97] T. Sjostrand, S. Mrenna and P. Z. Skands, JHEP 0605, 026 (2006) [hep-ph/0603175]. 

[98] J. Conway, PGS - Pretty Good Simulation of high energy collisions, http://physics.ucdavis.edu/ 
conway /research/software / pgs / pgs4-general.htm. 

[99] LHC SUSY Cross Section Working Group, https://twiki.ccrn.ch/twiki/bin/view/LHCPhysics/SUSYCrossSections 
[100] W. Beenakker, M. Klasen, M. Kramer, T. Plehn, M. Spira and P. M. Zcrwas, Phys. Rev. Lett. 83, 

3780 (1999) [Erratum-ibid. 100, 029901 (2008)] [hcp-ph/9906298]. 
[101] LHC Higgs Cross Section Working Group, S. Dittmaier, C. Mariotti, G. Pas- 
sarino, and R. Tanaka (Eds.), Handbook of LHC Higgs Cross Sections: 1. Inclu- 
sive Ohscrvables, CERN-2011-002 (CERN, Geneva, 2011), arXiv: 1101 . 0593 [hep-ph]. 



41 



https: / / twiki.cern.ch / twiki/bin /view/LHCPhysics/CrossSections 
[102] T. Junk, Nucl. Instrum. Meth. A 434, 435 (1999) [hep-ex/9902006]. 

[103] [ATLAS Collaboration], "Search for gluino-mediated scalar top and bottom quark production in final 
states with missing transverse energy and at least three 6-jets with the ATLAS Detector," ATLAS- 
CONF-2012-058 (2012). 

[104] [CMS Collaboration], "Search for supersymmetery in final states with missing transverse momentum 

and 0, 1, 2, or > 3 b jets in 8 TeV pp collisions" CMS-PAS-SUS-12-016 (2012). 
[105] S. Chatrchyan et al. [CMS Collaboration], arXiv: 1205.3933 [hep-ex]. 

[106] [CMS Collaboration], "Search for supersymmetry in events with same-sign dileptons" CMS-PAS-SUS- 
12-017 (2012). 

[107] S. Chatrchyan et al. [CMS Collaboration], JHEP 1206, 169 (2012) [arXiv:1204.5341 [hep-ex]]. 

[108] E. Contreras-Campana, N. Craig, R. Gray, C. Kilic, M. Park, S. Somalwar and S. Thomas, JHEP 

1204, 112 (2012) [arXiv:1112.2298 [hep-ph]]. 
[109] G. Aad et al. [ATLAS Collaboration], arXiv: 1204.6736 [hep-ex]. 

[110] [ATLAS Collaboration] Search for direct top squark pair production in final states with one isolated 
lepton, jets, and missing transverse momentum in sqrt(s) = 7 TeV pp collisions using 4.7 ifb of ATLAS 
data. ATLAS-CONF-2012-073 [2012] 

[111] [CMS Collaboration] Higgs to gamma gamma, Fermiophobic. CMS-PAS-HIG-12-022 [2012] 

[112] G. Aad et al. [ATLAS Collaboration], Phys. Rev. Lett. 108, 111803 (2012) [arXiv:1202.1414 [hep-ex]]. 

[113] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710, 403 (2012) [arXiv:1202.1487 [hep-ex]]. 

[114] [CMS Collaboration] Search for supersymmetry in events with photons and missing energy. CMS- 
PAS-SUS-12-018 

[115] [CMS Collaboration] "Search for a Higgs boson decaying into two photons in the CMS detector" 

CMS-PAS-HIG-11-010 (2011). 
[116] U. Baur, T. Plehn and D. L. Rainwater, Phys. Rev. D 69, 053004 (2004) [hcp-ph/0310056]. 
[117] R. Contino, M. Ghezzi, M. Moretti, G. Panico, F. Piccinini and A. Wulzer, arXiv: 1205.5444 [hep-ph]. 
[118] [CMS Collaboration], "Status of b-tagging and vertexing tools for 2011 data analysis" CMS-PAS- 

BTV-11-002 (2012). 

[119] [ATLAS Collaboration] Expected photon performance in the ATLAS experiment. ATL-PHYS-PUB- 
2011-007 (2011) 

[120] [CMS Collaboration] Top pair cross section in e/mu+jets at 8 TeV CMS-PAS-TOP-12-006 [2012] 
[121] [CMS Collaboration], "Search for the fermiophobic model Higgs boson decaying into two photons" 

CMS-PAS-HIG-12-002 (2012). 
[122] G. Aad et al. [ATLAS Collaboration], arXiv:1206.0756 [hep-ex]. 

[123] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710, 91 (2012) [arXiv:1202.1489 [hep-ex]]. 
[124] G. Aad et al. [ATLAS Collaboration], arXiv:1206.5971 [hep-ex]. 

[125] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 713, 68 (2012) [arXiv: 1202.4083 [hep-ex]]. 



42 



[126] G. Aad et al. [ATLAS Collaboration], arXiv:1207.0210 [hep-ex]. 

[127] S. Chatrchyan et al. [CMS Collaboration], Phys. Lett. B 710, 284 (2012) [arXiv:1202.4195 [hep-ex]]. 

[128] A. Arvanitaki and G. Villadoro, JHEP 1202, 144 (2012) [arXiv:1112.4835 [hcp-ph]]. 

[129] K. Blum, R. T. D'Agnolo and J. Fan, arXiv:1206.5303 [hep-ph]. 

[130] M. R. Buckley and D. Hooper, arXiv:1207.1445 [hep-ph]. 

[131] B. Fcigl, H. Rzehak and D. Zcppcnfeld, arXiv:1205.3468 [hep-ph]. 

[132] B. Fcigl, H. Rzehak and D. Zcppcnfeld, Eur. Phys. J. C 72, 1903 (2012) [arXiv: 1108.11 10 [hcp-ph]]. 
[133] M. Lisanti and N. Weincr, Phys. Rev. D 85, 115005 (2012) [arXiv:1112.4834 [hcp-ph]]. 
[134] D. Curtin, P. Jaiswal and P. Meade, arXiv:1206.6888 [hep-ph]. 

[135] [ATLAS Collaboration] Search for the Associated Higgs Boson Production in the WH -> WWW* -> 
Ivlvlv Decay Mode Using 4.7 ftrl of Data Collected with the ATLAS Detector at y/s = 7 TcV 
ATLAS-CONF-2012-078 (2012) 

[136] [ATLAS Collaboration] "Search for direct production of charginos and neutralinos in events with three 
leptons and missing transverse momentum in sqrt(s)=7 TeV pp collisions with the ATLAS detector", 
ATLAS-CONF-2012-077 (2012) 

[137] [CMS Collaboration] "Search for the standard model Higgs boson produced in association with W or 
Z bosons, and decaying to bottom quarks for ICHEP 2012" CMS-PAS-HIG-12-019 (2012) 

[138] T. Aaltoncn et al. [CDF Collaboration], Phys. Rev. Lett. 102, 221801 (2009) [arXiv:0903.2618 [hep- 
ex]]. 

[139] J. M. Buttcrworth, A. R. Davison, M. Rubin and G. P. Salam, Phys. Rev. Lett. 100, 242001 (2008) 
[arXiv:0802.2470 [hep-ph]]. 



